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Abstract 

We present a treatment of many-body Fermionic systems that facilitates an 
expression of the well-known quantities in a series expansion in Ti. The en- 
suing semiclassical result contains to a leading order of the response function 
the classical time correlation function of the observable followed by the Weyl- 
Wigner series, on top of these terms are the periodic-orbit correction terms. 
The treatment given here starts from linear response assumption of the many- 
body theory and in its connection with semiclassical theory, it makes no as- 
sumption of the integrability of classical dynamics underlying the one-body 

quantal system. Applications of the framework are also discussed. 
PACS numbers:05.30.-d 05.45. +b 
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1. Introduction 

Semiclassical framework for single-particle systems is presently in an advanced stage 
]1].|2|]. For both the integrable and chaotic dynamics, we now understand the semiclassical 
quantization. However, for intermediate and mixed dynamical scenario, we still lack a 
convincing theory. Another situation of a great practical relevance arises in the systems 
comprising of many bodies. That the spectral fluctuation characteristics of these systems 
are modelled in much the same way as for chaotic systems with lesser degrees of freedom has 
been shown quite recently @ . It is very important to note that complex systems do possess 
features, mainly associated with the generic nature of the thermodynamic limit, which are 
completely absent from the systems with fewer freedoms, friction being one of the examples. 
There are other systems which one would like to understand semiclassically. Some notable 
examples are of quantum dots 0, metallic clusters || and their optical properties 0; also, 
nuclear physics at high spins presents us with opportunity to test various semiclassical 
ideas in many-body theory. There has been an attempt based on the semiclassical limit of 
the time-dependent Hartree-Fock equation where one can establish a connection between 
the strength function given by the Vlasov equation and the corresponding quantum function 
in the limit of large quantum numbers. Analogous to the many-body treatment based on 
random-phase approximation (RPA), the linearized Vlasov equation gives an integral equa- 
tion for the particle-hole propagator in terms of classical propagator in the static mean field 
and two-particle Coulomb interaction. Thus one needs to evaluate the classical propagator 
in the mean field and then solve the integral equation. 

There is a lot of evidence that periodic orbits in the semiclassical trace formulae describe 
the shell effects in metallic clusters IJlQfl . Of course, these have been preceded by nuclear 
physicists |TI|] by more than a decade. In fact, the periodic-orbit theory ||12|| , which, in 
principle, should be rendered useless for practical purposes due to a huge number of technical 
problems associated with finding periodic orbits and summing a conditionally convergent 
series, gets encouragement from observations that only a few periodic orbits are enough in 
many situations 
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We present here a semiclassical treatment of the response function. Response function 
is essentially the imaginary part of generalized dynamical susceptibility which, in turn, is 
related to the correlation commutator. In the section 2, we present the general discussion 
from the linear response theory to arrive at a quantity for which semiclassical expression 
can be written. The key point of this section is to show that the most important quantity 
is a time correlation function. In the section 3, we present semiclassical treatment of the 
system perturbed by an external influence. This enables us to write the relevant results of 
sections 2 and 3 in a semiclassical expansion. We will see that the leading term is indeed the 
two-time correlation function averaged over the phase space, followed by the Weyl-Wigner 
series which has the periodic-orbit corrections. On our way, we take account of the Fermionic 
nature of the particles. 



2. Response Function and the Time Correlation Function 

We present here a discussion on the response of a system in the presence of an external 
field. We work with susceptibility |T3| and express it in terms of trace over single-particle 



states of two-time correlation function of the observable. This quantity is also known as the 
polarization propagator fd| . 

If a system described by the Hamiltonian H is externally disturbed by an field F ext (t), 
then the total Hamiltonian is 

H T = H- QF ext {t) (1) 

where Q is an observable, an example could be magnetization in the context of spin systems, 
or, an electric dipole operator in an example involving photoabsorption, and so on. The 
response function can be written as the imaginary part of the dynamical susceptibility, 

r(t,t') = (2h)- 1 ({Q(t),Q(t')]) 

= I ^e^f'H. (2) 

The angular brackets denote the expectation value and the square brackets denote the 
commutator. Setting an initial time to and the final time to t, we can write 



X "(t) = (2h)- 1 \(Q(t)Q(0)) - (Q(O)Q(t)) 



(3) 



where 



Q(t) = exp(iHt/h) Q exp(-iHt/h) 



(4) 



and where (•) denotes the average over the initial state of the system which is, for instance, 
the thermal state wherein 
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(5) 



with 



Z(P) =tr exp(-/3if), 



(6) 



if the system was in contact with a thermal reservoir at temperature T = l/(3k,B in the 
period preceding the interaction with the external field, F ext (t). If the system was in its 
pure ground state |$ ) then the average is over this many-body eigenstate of H, i.e., 



(•) = ($ | ' |$c 



(7) 



In this case, we can write x" (t) as 



x"(t) = {2%r 1 [<*o|Q(*)Q(0)|*o> - (%\Q(0)Q(t)\% 

= {2h)~ l [e iEot / h {%\Qe- iAt / n Q\%) - e - iE ° t/h (%\Qe tlit/n Q\<f> 
= (2Tiy 1 Y J [e tEot,n \{%\Q\^ n )\ 2 e- iEnt/h - e- iEot ' h \ (%\Q\$ n ) \ 2 e iE ^ h ] (8) 

n 

where {|$ n )} denote all the many-body eigenstates of the isolated-system Hamiltonian, 



H\$ n )=E n \$ n ), 

with E n > E . 

Upon Fourier transformation, 

X"(uj)=7ij2 \(®o\Q\$n)\ 2 [5(hu + E -E n ) -8(hu-E + E n )} 

n 

1 1 



(9) 



E l<*0|Q|*r. 



| 2 3 



hcu- E + E n + i0 + hu + E - E n + i0 + 



(10) 



where we used the identity 

— 1 =r(-) - t7c 5{x), (ii) 

V denoting the Cauchy principal value. Methods have been developed to evaluate semiclassi- 
cally such expressions as above |]J . However, we should recall that the system is many-body 
so that such methods would require the search for classical orbits of the many-body system. 
Simplification arises by taking account of the Fermionic character of the system which allows 
a reduction of the problem to one-body Hamiltonian in an effective potential determined 
by, for instance, the Hartree-Fock method. We shall make this simplifying assumption here 
and restrict our system to a set of uncoupled one-body Hamiltonians. A similar assump- 
tion is carried out on the coupling operator for which we assume the same form as for the 
Hamiltonian, 

TV 

# = E ^ 

i=l 
N 

Q = E Qi (!2) 

where hi and are one-body operators. Let us denote the one-body eigenstates as 

h\<f> a ) = e a \<t> a ). (13) 

When the Fermionic many-body system is in its ground state, all the one-body eigenstates 
are occupied up to the Fermi energy, so that the energy of the ground state is 

N 

E = E (14) 

£ a( f F 

In any many-body excited state, |$ n ), at least one 1-body level above the Fermi energy is 
occupied. We may denote such excited states by the list of the occupied 1-body levels, or 
equivalently, by the list of the one-body levels for which the occupation is different with 
respect to the ground state |$o)> having in mind that such states are antisymmetric for an 
exchange of two Fermions: 



i$ ) = iiiii-ii ; eF oo. ..o) 

|$ n ) = 1 111...1110 - F 0100...). (15) 

For operators which are sums of 1-body operators as assumed in (12), we obtain the result 
that the matrix elements (<& |Q|<l> n ) are non- vanishing only for the states $ n which differ 
from the ground state $ by one excitation. These states have a hole in state a (e a (e F ) and 
a particle in state b (e^ep). For these states the matrix elements are thus 

($ |Q|$n) = (<f> a \q\<f>b). (16) 

Moreover the energy of the excited state is 

E n = E -e a + e b (17) 
under the assumption that e a {ep{^b- After some standard manipulations, we can write 

\{^\Q\®n)\ 2 



E 



hu- E + E n + i0+ 



-g ft ,~;,:0+ [l-e(e r - et )]e(e F -0, (!8) 

where the Heaviside step function, ©(■) takes care of the aforementioned restriction on the 
location of the 1-body states <fi a and (fib with respect to ep. 

In the presence of a thermal reservoir at temperature T[— l/(/3k B )] described with 
canonical density matrix, we get that 

(Q(t)Q(0)> = | tr e-^e^Qe-^Q 

Zj 

= '£e X p[iE m (t + iph)/h] exp(-tE n t/h) | ($ m \Q\Q n ) | 2 . (19) 

mn 

Thus the dynamical susceptibility is 

Similar assumptions as before enable us to reduce this expression to the 1-body system. 
Now the probability to find a state |$ m ) in which \<fi a ) is occupied is given by the Fermi- 
Dirac distribution at energy e a . But is related to |$ m ) by the fact that \<fib) must be 
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unoccupied (see above). So, we have a joint probability that \<p a ) is occupied and \4>b) is 
unoccupied, hence the propagator is 

FD _ FD 

n^-ZKt.m*?*^^ (2D 

where 

p FD = (22) 

denotes the Fermi-Dirac probability or mean occupation number, /i is the chemical potential 
fixed by the total number of Fermions in the system. Eq. (21) gives therefore the Fourier 
transform of the two-time correlation function. 

At zero temperature, the Fermi-Dirac distribution becomes a Heaviside step function, so 
the expression (21) becomes 

ab nuj — e a + e b + 2ir 

which is identical with the expression previously derived. In the foregoing discussion of this 
section, we have started from the time correlation function (actually the commutator) and 
written an expression for the response. Now that the susceptibility has been reduced to a 
1-body expression, we can treat the 1-body system semiclassically in terms of periodic orbits 
in the effective Hartree-Fock Hamiltonian h. 

It is well to recall that x" (^) has a physical interpretation in terms of energy dissipa- 
tion which is related to the fluctuations in the frequency spectrum obtained as the Fourier 



transform of the time correlation function. Indeed the fluctuation-dissipation theorem [13 
lies just in realizing that connection. 

Employing the identity (11) in (21), and expanding the difference between p FD 's in the 
resulting expression, we obtain 

X" (lu)=itJ2\ (<l>a\q\(t>b)\ 2 [p FD (e a ) - p FD (e a + huo)} 5(huj + e a - e b ) 

a,b 

00 (flLj) n f)n n FD 

= " 7r E^pE ~Lr^) mWb)? S(hco + e a - e h ) 

n=l ' a,b 



where we introduced the expression 

QripFD 



/» = 27rft]T -J^-(e a ) |(0 a |g|0 b )| 2 5(^ + e a -e 6 ) = J dt f n (t), (25) 

a,b 

which is the Fourier transform of the time correlation functions 

an FD 

f n (t) = tr -^-(h) e ah ' h q e~ ith ' n q. (26) 

Moreover, we can rewrite these time correlation functions in terms of a single time correlation 
function as follows 

/an FD 
dz^riz) W), (27) 

with 

C e (t) = tr S(e - h) q(t) q(0) and q(t) = e ith/h q e~ ltk/n . (28) 

In summary, we have so far reduced the many-body dynamical susceptibility to an expression 
involving the time correlation function C t (t) of the 1-body effective dynamics: 

We will see in the next section that this expression is particularly suited for developing 
semiclassical expansions. 

3. Semiclassical Expressions of Response 

According to the preceding section, we need to write the semiclassical expression for the 
time correlation function C e (t) in order to obtain semiclassically the response function. To 
this end, we consider the correlation function in the form: 

C e (t) = tr 6{e - h)X = tr 6{e -h)A + i tr 5(e - h)B, (30) 

where X = X(t) = q(t)q(0) is a 1-body operator. This operator is nonHermitian but can 
be decomposed as X = A + iB in terms of two Hermitian operators: A = (X + X^)/2 and 
B = (X- A"t)/(2i). 



Methods have been obtained to evaluate semiclassically such expressions in terms of 
periodic orbits fllSfl. To motivate this method, we begin by observing [|1^] that such an 
expression involves the matrix elements of A (and B) over the eigenstates of h. These matrix 
elements can be obtained at the first-order perturbation theory of a perturbed Hamiltonian 
h(X) = h + XA. Assuming the eigenvalue problem for h(X) to be solved, the matrix elements 
of A may thus be obtained in terms of derivatives of eigenvalues of the perturbed Hamilto- 
nian, e n (A), with respect to A. Now the periodic-orbit theory by Gutzwiller can be used to 
calculate the diagonal matrix elements. 

Each term of the correlation function C e (t) can therefore be expressed as 



tr S(e - h)A 



xs tr 



.4 



7T 



e - h + i0 + 



Id - 

- 3 tr — log(e - h - XA + i0~ 

7T oX 



A=0 



(31) 



Comparing with the following indentity 
1 d 



7i de 



9 tr log(e - h — XA + i0 + ) = tr 5(e - h - XA) 

= |iV(e;A), 



(32) 



which is the derivative of the staircase function iV(eA) with respect to the energy, we arrive 
at the relation, 

dN{e; A) 



tr S(e - h)A 



OX 



A=0' 



(33) 



which gives each term of the correlation function as the derivative of the staircase function 
of the perturbed Hamiltonian h(X) with respect to A. 

On the other hand, one has the well-known semiclassical expression: 



N(e;X) 



(2nh) 



&[e-h w (X)} + o(h-f +1 ) 



00 1 sin 

+ EE 

p r- 



jrS p (e; A) - rf u p 



+ o(h), 



(34) 



^rir |det[m;(A) - l]) 1 / 2 
where m p (A) is the monodromy matrix governing the stability of the classical periodic tra- 
jectory, p; u p is the Maslov index of the trajectory and hw is the Weyl-Wigner transform 
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of the Hamiltonian h(X) = h + XA. We notice that the above periodic orbits p are those 
of the perturbed Hamiltonian. We assume here that the periodic orbits of hw{X) deform 
continuously to the periodic orbits of hw(X = 0) for A small enough. For the purpose of 
differentiating this expression with respect to A, we use the following classical formula which 



gives the derivative of the action of the periodic orbits |19 

dh w (X) 



dt 



dt A 



w • 



(35) 



OX J OX 

where the integrals go around the periodic orbit, p, and Ay/ is the Weyl-Wigner transform of 
the operator A. Taking the derivative with respect to A and adding both terms composing 
the correlation function, we finally obtain 



tr 5{e - h)X 



d xd p 

(2^)7 Xw 6{e ~ kd) 



0{h 



I 00 cos 

+ —YY — 

nh^f^ |det(mr 



0| 1/2 



dt X> 



w 



0{h°) 



(36) 



p r=l I -**">\—p 

where h c i = hw(X = 0) and the periodic orbits are those of h c i. 
For the case like ours, when 



X = exp(-ht) q exp(--ht) q 
a k 



the Weyl-Wigner tranform can be written as 



(37) 



X w (x, p) 



. e 2 qw 



(e* ht ) w e^Hq w e^ k [(e^ ht ) w 
exp ^— C c itj qw qw + (Weyl corrections 
qw(t)qw{0) + (Weyl corrections) 



(38) 



where C c i = {h, ■} is the classical evolution operator, and the operator A is standard [g0| in 
the Weyl-Wigner expansions defined as 



A 



d d d d 



(39) 



dp <9x <9x dp 

and qwif) = Qw [^'( x ? p)] is the value of qw at the current point $*(x, p) of the trajectory 



from the initial condition (x, p) of the Hamiltonian flow $* of h 



cl- 
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Putting all the pieces together, we can now write the semiclassical expansion for the 
correlation function as 

C e (t) = tr 5(e - h) q(t) q(0) 

(2n%)f 5 ^ e ~ P ^ qw<y *' P ^ ( e ~ £ci< ^)( x ' P) + ( We ^ corrections) 
1 cos (tS„ — r^uS) r 

+ UK 5 idetK-oV I dT + () + 0(h ) - (40) 

We now go back to Eq. (29) where we see that we still need to evaluate the Fourier transform 
of the correlation function C e (t). 

Let us evaluate the first term with n — 1 of the susceptibility. This should give the 
dominant term of the Weyl series in powers of the Planck constant H as well as a sum 
over periodic orbits with the smallest power in h which gives therefore the most important 
contribution. The other periodic-orbit sums would involve smaller amplitudes with higher 
powers in h. Replacing with the previous semiclassical result, we get 

Q p FD 



» = ~Jde *L-(e) J dt C e (t) + 0(h) 

= -^j dt J d e p FD [h d (x,p)} q w (x,p) [e^-^q w ] (x, p) + 0{h~^) 

r ^ cos S v — r^u v ) r r 



2nh 

Assuming that q = q* is a Hermitian operator, its Weyl-Wigner transform is a real function 
of (x, p). Along a periodic orbit, we have the Fourier series 

^ A 27m \ 

Qw{r) = }^ %,n exp U^— r , (42) 

n=-oo V 1 P J 

where T p is the classical period and q p - n = Qp >n - We obtain that 

27TO 

Except in special systems like the harmonic oscillators, the periods T p (e) vary with the 
energy e. Because of the presence of the Dirac distribution, we can perform the last integral 
over e by using 
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dte™ fdrq w (T)q w (r + t)=27rT p (e) £ |g P ,„(e)| 2 6 

J P n=-oo 



(43) 



UJ — 



2im 



T 2 



2irnd e T p 



(44) 



where e = e Ptn (uj) is the energy obtained by solving the implicit equation T p (e) = 2nn/uj. 
We here assumed that the periods vary monotonously with energy, otherwise we find a su 
of terms instead of a single term in the right-hand member of the previous formula. Noting 
that 



dS pn d . . 



rp3 
P 



2nnd t T p ' 



(45) 



we have that 



r r d^ 

J dt jdr qw(r) q W (r + t) = 2tt]T |g p , n ( e )| 2 \^\ 5 [e - e p , n (uj)] , 



(46) 

where S p>n = S p [e p<n (uj)]. 

Substituting in the previous relation, we finally obtain the absorptive part of the dynam- 
ical susceptibility as: 

d*xdfp 



UJ 



d,t 



UJ 

h 



COS (jfiS Pjn '"'f^p) 



e (iu,-C cl )t qw 



p,r,n 





dS Pt n 


\p,n \q.p,n\ 


duo 



+ o(h°). 



x, P ) + o(h- f+1 ) 

(47) 



det( m ; jn -l)|V2 

where all the expressions with indices p, n are functions of the frequency u> through their 
dependency on the energy e = e Pin (u). Eq. (47) is the central result of this work. 

To interpret this result, let us go back to the time correlation function C € (t) which is 

tr 5(e - h) exp(iht/h) q exp(-iht/h) q = ^- f dr e ter/h tr e ih(t - r)/h q e~ lkt/h q. (48) 

2n J 

In this expression, the trace can be rewritten in terms of the Green propagators as 

tr exp[ih(t — r)/h] q exp[—iht/h] q 
= J dxdx' G(x',x;T-t) g(x) G(x,x';t) g(x') 
= ^2 J dxdx.' A c (x', x; r — t) exp[iW c (x', x; t — t)/h] q(x) 

c,c' 

xA c /(x,x';t) exp[iW c ,(x,x.';t)/h] q(x.') (49) 
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where in the last step, we have written down symbolically the semiclassical expressions for 
the Green propagators. We have here assumed that the operator q depends only on the 
position operator, i.e., q = <?(x). The sum in the above expressions is over classical paths, c 
and d. Notice that these are two different paths in general as they arise from two different 
propagators. The expression interprets for itself - from the point x to x', we propagate along 
the path d for time t, and then we propagate back to the point x along c in time r — t. 
Precisely, 



L being the Lagrangian of the system. The total phase in the correlation function is the 
sum W c + W c i of the two actions. By stationary phase method, it can be shown that the 
two segments of the paths c and d, in fact, form a periodic orbit of period r. This is the 
mechanism by which periodic orbit corrections appear in the main result above. 

Classical contribution comes from the paths of zero length, i. e., r = 0, and this is 
precisely the leading term in our expression, also called the quasiclassical term. Thus, 
in any calculation of interest, we need to find the classical correlation function, the Weyl 
corrections, and then the periodic-orbit corrections. If the energy spectrum of the system is 
discrete, we notice that the dynamical susceptibility (24) is a sum of Dirac peaks centered 
on each of the Bohr frequencies 00 = (E m — E n )/%. In the semiclassical formula (47) the 
first term due to the paths of zero length contribute to the dynamical susceptibility by a 
continuous function of frequency u which may be interpreted as the average or background 
of the discrete quantum susceptibility. The second series of terms of (47) are due to the 
periodic orbits which contribute by oscillatory functions of the frequency uo. Their sum is 
supposed to reproduce the Dirac peaks of the discrete quantum susceptibility 

We notice that, in contrast with the well-known Gutzwiller formula for the level density, 
the quasiclassical term here involves a classical time correlation function. As a consequence, 
the spectral properties of the classical Liouvillian t c i will intervene in the evaluation of this 




(50) 
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term (as well as of the following Weyl corrections). In this regard, we mention that a periodic- 
orbit theory of the classical Liouvillian has been developed p!7JI^| in terms of classical 



zeta functions which resemble but differ from the Selberg-Gutzwiller quatum zeta function. 
Consequently, we have here an example where the semiclassical evaluation of a quantum 
property, namely involves the spectral properties of the classical Liouvillian along 

with the Gutzwiller quantum periodic-orbit sum. It is remarkable that both these classical 
and semiclassical terms coexist within the same formula, which points to the result that the 
classical periodic-orbit theory is only one part of the underlying quantum dynamics. We 
believe that this result is general in the semiclassical evaluation of dynamical susceptibilities, 
which include in particular the ac conductivity. 

It turns out that the spectral properties of the Liouvillian differ considerably whether 
the classical system h c i is integrable or chaotic. 

For integrable systems, the Liouvillian spectrum is formed by discrete classical frequen- 
cies, Cdipm = ifi m ^ m , on each ergodic component which are the invariant tori. The expres- 
sion would then be similar to the one obtained in J7J. The periodic-orbit corrections should 
be given by the Berry- Tabor semiclassical theory under such circumstances . 



For chaotic systems, the classical time correlation function in the quasiclassical term 
decays if the system is mixing on the energy shell. Recent works have shown that this decay 
may be developed in terms of the complex singularities of an analogue of the Liouvillian 
resolvent |T7| , p!8|p3| . When the complex singularities are poles they are called Pollicott-Ruelle 
resonances and are associated with exponential decays. If the spectral decomposition of the 
Liouvillian dynamics in terms of the Pollicott-Ruelle resonances is known the quasiclassical 
term can be further reduced as a sum over the Pollicott-Ruelle resonances over each energy 



shell P"7|,|i8ira 



4. Concluding Remarks 

We have described a systematic procedure which lets us write important quantities in 
many-body theory in a semiclassical expansion. All our considerations are restricted to 
many-body systems reducible to uncoupled one-body systems, we believe that it is an inter- 
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esting first step. We have shown the central role played by the time correlation functions. 
On the one hand, they are related by a Fourier transform to the response function; on the 
other hand, they facilitate useful expansions in terms of h. The fact that the Fermionic 
nature of the particles has been taken into account is an important aspect. The other rea- 
son for an interest in the above considerations stems from the fact that no assumptions are 
made about integrability of the underlying classical one-body dynamics. It may be recalled 
that in a recent attempt along similar lines @, complete integrability was assumed. More- 
over, in their analysis, the Berry- Tabor formula for periodic orbit corrections has been used 
which fails for the interesting case of the harmonic oscillators. Since we now know that 
the Gutzwiller trace formula is exact for harmonic oscillators |J, and is also applicable for 
chaotic systems, we believe that our analysis extends to more general systems. 

There are many applications of such a formalism. In any discussion of relaxation phe- 
nomena, time correlation functions play an important role. To understand photoabsorption 
cross-section of atoms and molecules, semiclassical treatment along these lines has been used 
PIPT|] . In metallic clusters and quantum dots, there is a growing interest in the semiclassical 
treatment of the absorption spectra and cross-sections. The understanding of the plasmon 
modes in different metallic clusters, their splitting etc. are believed to have a semiclassical 
interpretation. However, for a systematic treatment of the collective modes, we need to 
consider two-body interactions, thus a Hamiltonian, 

H = y £h i + X y £v ij (51) 

i ij 

where denotes the two-body term. What we need to do here is a perturbative treatment 
of the two-body interaction by assuming that A 1 which, in turn, would give rise to 
collective modes. 

In a random-matrix framework, one can also study the correlation functions, and it 
has been recently shown that the time correlations depend upon the co-dimension of level 



crossing |24j|. A comprehensive understanding of all the inter-relations will be one of the 



aims we look forward to. 
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